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Abstract 

The role of physical barriers in promoting population divergence and genetic struc- 
turing is well known. While it is well established that animals can show genetic 
structuring at small spatial scales, less well-resolved is how the timing of the ap- 
pearance of barriers affects population structure. This study uses the Panama Canal 
watershed as a test of the effects of old and recent riverine barriers in creating pop- 
ulation structure in Saguinus geoffroyi, a small cooperatively breeding Neotropical 
primate. Mitochondrial sequences and microsatellite genotypes from three sam- 
pling localities revealed genetic structure across the Chagres River and the Panama 
Canal, suggesting that both waterways act as barriers to gene flow. F -statistics and 
exact tests of population differentiation suggest population structure on either side 
of both riverine barriers. Genetic differentiation across the Canal, however, was 
less than observed across the Chagres. Accordingly, Bayesian clustering algorithms 
detected between two and three populations, with localities across the older Chagres 
River always assigned as distinct populations. While conclusions represent a prelim- 
inary assessment of genetic structure of S. geoffroyi, this study adds to the evidence 
indicating that riverine barriers create genetic structure across a wide variety of taxa 
in the Panama Canal watershed and highlights the potential of this study area for 
discerning modern from historical influences on observed patterns of population 
genetic structure. 



Introduction 

The distribution of genetic variability across geography is 
affected by multiple biotic and abiotic factors (Loveless and 
Hamrick 1984; Avise 2004), including mode of reproduction, 
vagility, philopatry, and geography. The role of physical bar- 
riers in promoting population divergence and structure is 
well known (Avise and Felley 1979; Preziosi and Fairbairn 
1992). While it is well established that animals can show ge- 
netic structuring at very small spatial scales (Selander 1970), 
the timing of the origin of those patterns has been more dif- 
ficult to discern because they may be the combined result 
of contemporary processes (Zellmer and Knowles 2009) as 
well as longer term historical events (Bowen and Avise 1990). 
Thus, the timing and appearance of physical barriers and how 
quickly these affect genetic structure in populations remains 
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a topic of interest (Matocq et al. 2000; Vandergast et al. 2007; 
Zellmer and Knowles 2009). 

Landscape features that constitute barriers vary among 
species. For instance, differences in elevation contribute to 
population differentiation in two amphibians, the blotched 
tiger salamander (Spear et al. 2005) and the Columbia spot- 
ted frogs (Funk et al. 2005). Conversely, pacific jumping 
mice {Zapus trinotatus) readily bound large topographic 
barriers, with gene flow explained more appropriately by 
habitat connectivity (Vignieri 2005). Even relatively new 
and small barriers can affect the population structure of 
animals. Epps et al. (2005) showed that recent (~40 yrs) 
anthropogenic barriers have caused a marked decline in 
genetic diversity, in a large vagile mammal, the desert 
bighorn sheep {Ovis canadensis nelsoni). Similarly, anthro- 
pogenic barriers that red grouse could theoretically cross 
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in one flight acted as barriers to gene flow (Piertney et al. 
1998). 

Bodies of water can promote genetic differentiation in a 
variety of terrestrial species. Sea lochs explain most of the 
genetic differentiation among populations of red deer in the 
Scottish highlands (Perez-Espona et al. 2008). Quemere et al. 
(2010) found that the Manankolana River was the major 
barrier to gene flow for the golden crown sifaka (Propithecus 
tattersalli) in the Daraina region of Madagascar. Bodies of 
water may also be barriers for volant animals, as reported 
by Meyer et al. (2009) for some bat species. Rivers can serve 
as barriers for amphibians, as reported for the alpine stream 
frog (Scutiger spp.) in the Hengduan Mountains of China (Li 
et al. 2009). 

Although riverine barriers have been implicated as barri- 
ers to gene flow in a variety of species, less is known about 
how the timing of the appearance of these barriers affects 
genetic structure. The Panama Canal is one of the largest 
modifications of the hydrographic landscape undertaken by 
humans. Because the addition of this major riverine bar- 
rier is well documented historically and the Panama Canal is 
embedded in the center of a dynamic watershed with older 
riverine barriers, it presents an ideal opportunity to test the 
influence of the timing of physical barriers in population ge- 
netic differentiation. Before the construction of the Canal, 
the Chagres River flowed along the Atlantic slope of the Isth- 
mus, while only small costal streams drained the Pacific slope 
(Meek and Hildebrand 1916, Fig. 2). Geologic studies indi- 
cate that the Panamanian Isthmus was formed by sometime 
in the Pliocene (Coates et al. 2004; Kirby et al. 2008), with 
little tectonic activity after the late Pliocene (Coates et al. 
2004) indicating that by then the topographic composition 
of the Panama Canal watershed was probably similar to that 
observed today. Thus, while the Chagres River had been a 
major part of the central isthmian basin in paleontological 
time, the completion of the Panama Canal in 1914 created a 
novel riverine barrier, which is expected to have affected the 
movements and, hence, gene flow in a multitude of species. 

This study uses the Panama Canal watershed to test the 
effects of an old and a recent riverine barrier in creating pop- 
ulation structure in Geoffroy's tamarin Saguinus geoffroyi, a 
small Neotropical primate (Fig. 1). Tamarins represent an 
appropriate study species for this study for several reasons, 
first, the Chagres River (the old river) was in place as a ma- 
jor riverine barrier well before the estimated divergence time 
(0.7 Mya, Evans et al. 1998) of Geoffroy's tamarin from its 
sister species (S. oedipus). Second, previous work has shown 
that riverine barriers are important for structuring primate 
populations in general (Wallace 1852), and for the diversifica- 
tion of tamarins in particular. Hershkovitz (1977) originally 
proposed that subspecies of S. fuscicollis were delineated by 
rivers based on pelage color variation, a hypothesis later sup- 
ported by analyses of craniofacial variation (Cheverud and 
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Figure 1. Geoffroy's tamarin in the Gamboa forest. Photo by Anand 
Varma. 

Moore 1990). Peres et al. (1996) showed that cytochrome b 
haplotypes corresponded to phenotypically distinct morphs 
of subspecies of S. fuscicollis on opposite sides of the Rio Jurua 
in Amazonian Brazil. Furthermore, the authors showed that 
gene flow (and intergradation of color morphs) increased to- 
ward the narrow headwater streams of the river, as predicted 
by Wallace (1852). Third, among the primates inhabiting ar- 
eas close to the Panama Canal (howler monkeys: Alouatta 
palliata and capuchin monkeys: Cebus capuchinus), tamarins 
may be most likely to exhibit rapid population differenti- 
ation in response to landscape changes. Tamarins exhibit 
high intragroup relatedness (Huck et al. 2005), suggesting a 
low dispersal rate and a greater likelihood of showing genetic 
structure. In contrast, howler monkeys inhabiting an isolated 
island in the Panama Canal watershed for > 90 years showed 
no genetic evidence of a population bottleneck (Milton et al. 
2009), suggesting a near panmictic mating pattern. 

To examine the role of rivers in creating population ge- 
netic structure at a small spatial scale, I sampled three popu- 
lations distributed across two prominent riverine barriers — 
the Panama Canal and the Chagres River — to test whether 
the age of a physical barrier to gene flow has an effect on the 
level population genetic structure. I predict that there will 
be significant differentiation between populations separated 
by the Chagres, whereas differentiation across the Panama 
Canal will be more modest, owing to its novelty as a bar- 
rier: whereas the individuals included in this study repre- 
sent approximately 20-40 generations since the construction 
of the Canal, the number of generations since the appear- 
ance of the Chagres is at least three orders magnitude higher. 
To test these predictions, I: (1) examine differences in ge- 
netic variability among sampling localities at mtDNA and 
microsatellite loci, (2) use f -statistics and analysis of molec- 
ular variance (AMOVA) at mtDNA and microsatellite loci 
to investigate population genetic structure between sampling 
localities, and (3) use Bayesian clustering algorithms to de- 
termine the number of likely populations. 
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Figure 2. Sampling localities in the Panama Canal watershed. The top panel shows the current watershed, after the construction of the Canal. Inset 
shows the haplotype parsimony network generated by TCS, with haplotypes color coded to correspond with sampling localities and proportionally 
sized according to haplotype frequency. The bottom panel shows the watershed before the construction of the Canal. 



Materials and methods 

Study sites 

I targeted three areas within the Panama Canal watershed 
(provinces of Panama and Colon, Republic of Panama) 



for sampling (Fig. 2): (a) The Soberania field site located 
within the boundaries of Soberania National Park, close to 
Camino de Plantation (9.076°, -79.659°). (b) The Gamboa 
field site, just outside the park boundaries, located in and 
around the rural town of Gamboa, Colon Province (9.118°, 
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-79.698°). (c) The Panama West locality (8.957°, -79.668) 
west of Panama City across the Panama Canal. I selected 
the Soberania and Gamboa field sites to provide a compar- 
ison across the Chagres River; these samples were collected 
as part of a separate study on cooperative breeding in S. 
geoffroyi (Diaz-Munoz 2011). I selected the Panama West 
locality to capture the potential effect of the Panama Canal 
as a barrier. I collected samples from Soberania National 
Park and Gamboa from field captures, whereas the samples 
from western Panama Province were obtained from museum 
skins collected by Dawson (1976) and housed at the Michi- 
gan State University Museum (Table 1). The use of museum 
samples imposed limitations on the study, specifically they 
represented a sample from a wider geographic area as com- 
pared to the other localities and represented a time period 
approximately 30 years earlier than field-collected samples 
(field samples: 2005-2009; museum samples: 1973-1974). 
However, I opted to use museum samples due to the fact 
that some areas directly across the Panama Canal from 
Soberania have limited accessibility (due to deposits of unex- 
ploded ordnance from military activities) and tamarin pop- 
ulations at some localities in Western Panama have been 
extirpated due to growing urbanization. Utilizing museum 
resources allowed a broader sampling than would have been 
possible with field efforts and allowed looking at both of 
the riverine barriers of interest, albeit in a preliminary 
way. 

Captures and sample collection 

I captured tamarins from Gamboa and Soberania using 
hand-activated live traps baited with bananas as described by 
Garber et al. (1993) or by blow-darting (Bioject, Blowguns 
Northwest, Richland, WA) with tranquilizer darts (Pnueu- 
Dart, Williamsport, PA). Individuals from the Soberania 
population were only captured in traps and were not anes- 
thetized. To prevent excessive stress, I limited handling time 
< 15 min and manipulations were constrained to mark- 
ing and sampling hair. Soberania individuals were released 
immediately at the capture site. I captured Gamboa adults 
using blow-darting; infants and juveniles were trapped be- 
cause they were too small to safely dart. In both trapping and 
darting, I anesthetized Gamboa individuals using Ketamine 
(7.5 mg/kg) and Zoletil (3.75 mg/kg Vibrac SA, Carros Cedex, 
France). Gamboa animals were anesthetized to enable collec- 
tion of morphological data for a separate study (Diaz-Munoz 
2011). I handled Gamboa individuals for 48 ± 14 min and 
placed them in a pet kennel for 3.67 + 2.13 h until fully re- 
covered. I measured respiration, heart rate, and internal body 
temperature of anesthetized individuals throughout handling 
procedures to monitor animal condition. To minimize po- 
tential injuries, I darted individuals at feeding stations that 
were eye-level above the ground. After darting, individuals 
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Table 1. Michigan State University Museum specimens sampled for this 
study. 



Catalog 
no. 


Sex 


Year 

collected 


Specific locality 


Latitude 


Longitude 


MR.22872 


M 


1973 


Cerro Cama 


9.01667 


-79.90000 


MR.22874 


M 


1973 


Vicinity of La 


8.88333 


-79.78333 








Chorrera 






MR.22878 


F 


1973 


1 5 km W of 


8.95000 


-79.70306 








Balboa 






MR.22875 


M 


1973 


1 S km W of 


8.95000 


-79.70306 








Balboa 






MR.22947 


M 


1973 


6.5 km NW of 


8.98197 


-79.61785 








Balboa 






MR.22885 


M 


1973 


5 km NE of 


8.96598 


-79.63393 








Arraijan 






MR.22891 


M 


1973 


3.3 km NE of 


8.97110 


-79.62878 








Arraijan 






MR.22889 


F 


1973 


2.5 km NEof 


8.95000 


-79.58635 








Arraijan 






MR.22998 


F 


1974 


3 km W of 


8.95000 


-79.56816 








Balboa 






MR.22989 


M 


1974 


8.5 km Wof 


8.95000 


-79.64850 








Balboa 






MR.22935 


F 


1973 


6 km SW of 


8.92924 


-79.61707 








Balboa 






MR.22994 


F 


1974 


8.5kmWSW of 


8.92762 


-79.58917 








Balboa 






MR.22923 


M 


1973 


3.5 km SW of 


8.89885 


-79.61811 








Balboa 






MR.22963 


F 


1973 


4 km ESE of 


8.91164 


-79.60525 








Arraijan 






MR.22907 


M 


1973 


9 km W of 


8.99156 


-79.60846 








Balboa 






MR.22902 


F 


1973 


9 km E of 


8.95000 


-79.61363 








Arraijan 






MR.22949 


F 


1973 


4 km E of 


8.93616 


-79.61640 








Arraijan 






MR.22934 


M 


1973 


8 km SW of 


8.95865 


-79.62900 








Ddl UUd 






MR.22895 


F 


1973 


7 km E of 


8.96730 


-79.60800 








Arraijan 






MR.22985 


F 


1973 


5 km ENE of 


8.95000 


-79.64396 








Arraijan 






MR.22980 


F 


1973 


2.5 km ENE of 


8.92059 


-79.63807 








Arraijan 






MR.2291S 


F 


1973 


6 km WSW of 


8.95000 


-79.59395 








Balboa 







were followed by two field assistants with a mesh net to catch 
anesthetized individuals that strayed from the feeding sta- 
tion. All capture and handling procedures were approved by 
the UC Berkeley Institutional Animal Care and Use Commit- 
tee and followed the guidelines of the American Society of 
Mammalogists (Gannon and Sikes 2007). 

Field collected samples included: (a) Hair samples plucked 
from the base of the tail and saved in coin envelopes and 
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stored dry and (b) Ear tissue collected from the pinnae using 
surgical scissors and stored in KNAlater and frozen at -20° C 
until extraction. Soberania animals were represented by hair 
samples and Gamboa individuals were represented by hair 
and tissue samples. To verify the reliability of microsatel- 
lite genotypes from hair samples, I genotyped all Gamboa 
individuals using tissue samples to corroborate the results 
obtained from hair samples. 

Museum skins were preserved as dry flat skins. I used sur- 
gical scissors to extract ca. 1 mm 2 piece of tissue from the edge 
of the flat skin. I stored tissue samples in an empty micro- 
centrifuge tube until extraction, which occurred within the 
week. A list of sampled individuals is included in Table 1. 

DNA extraction 

I extracted Genomic DNA using Qiagen DNA Micro kits (Qi- 
agen, Valencia, CA), according to manufacturer instructions 
for each sample type. I soaked museum tissue samples in 70% 
ethanol for 24 h prior to extraction. I quantified DNA yield 
using a NanoDrop Spectrophotometer (Thermo Fisher Sci- 
entific, Waltham, MA). I extracted DNA from hair and mu- 
seum skin samples in a "clean" room dedicated to low-copy 
sample extractions. Both sampling and extraction negative 
controls were used to monitor for possible contamination at 
every step of the genetic workflow. 

Mitochondrial sequencing 

Mitochondrial sequences are extensively used in population 
genetic studies (Avise 2004) owing to their uniparental in- 
heritance and their relatively rapid sequence evolution, es- 
pecially in the mitochondrial control region (Hoelzel et al. 
1991). I amplified a 1080-bp fragment of the mitochon- 
drial control region using primers designed to be genus- 
specific for Saguinus (Table 2; Cropp et al. 1999). To design 
primers, which would amplify short fragments suitable for 
skin and hair samples, I used control region sequences de- 



Table 2. Mitochondrial DNA primers used in study. 



Primer 
Name 


S'-3' sequence 


Reference 




SCJB 


TTG GTTATGTAATTAGTG C 


Cropp et al. 


1999 


SCJ1 


GAGCGAGAATACTAGTAGAAG 


Cropp et al. 


1999 


464 


TGAATTGGAGGACAACCAGT 


Cropp et al. 


1999 


SCJ4 


GCACTAATTACATAACCAA 


Cropp et al. 


1999 


282 


AAGGCTAGGACCAAACCT 


Cropp et al. 


1999 


SCJ2 


ACCCTTCAGAGAATAAACTTA 


Cropp et al. 


1999 


SGDL1-F 


GCACACGACTACCAAGCAAGATTATGA 


This study 




SGDL1-R 


GGGTGGGGTGGGGACCAAGA 


This study 




SGDL6-F 


TCATCAGCATTGCCGGGAGGC 


This study 




SGDL6-R 


TGGTAGGCTAGGGGTATGTGGGG 


This study 




SGDL7-F 


ACCCAAAAATCACCACCCTAAGCG 


This study 




SGDL7-R 


TGGGGTTGTGACTGCCCATCT 


This study 





rived from tissue-extracted DNA from the Gamboa popula- 
tion. I amplified and sequenced the majority of the control 
region (ca. 1300 bp) — using "universal" primers MVZ 121/70 
and MVZ 123/106 derived from other studies (Kocher et al. 
1989; Palumbi 1996)— to design primers SGDL1-F, SGDL1- 
R, SGDL6-F, SGDL6-R, SGDL7-F, SGDL7-R using Primer3 
(Rozen and Skaletsky 2000), as implemented in Geneious 
4.8.5 (Biomatters, Auckland, NZ). I performed PCR reac- 
tions on ABI 2720 (Applied Biosystems, Foster City, CA) or 
BioRad iCycler (Bio-Rad, Hercules, CA) thermocyclers using 
fluorescently labeled primers. Cycling conditions were 94°C 
for 4 min, 94°C for 1 min, 52°C for 1 min, 72°C for 75 sec, 
repeated 30 times; 72°C for 10 min. Polymerase chain reac- 
tion volume was 10 yu.1 with 40 ng of genomic DNA, 1 fil 
of 10X PCR Buffer (Applied Biosystems), 2.5 mM MgCl2, 
0.8 ill 10 mg/mL BSA, 0.4 mM of each DNTP, 3 pM of each 
primer, and 0.5 U of Taq polymerase (Invitrogen, Carlsbad, 
CA). I confirmed amplification via TBE (Tris-borate-EDTA) 
agarose gel electrophoresis and cleaned products of ampli- 
fication using Exo-SAP IT (Affymetrix, Cleveland, OH). I 
fluorescently labeled PCR products utilizing ABI Big Dye 3.1 
(Applied Biosystems) and sequenced amplicons in an ABI 
3730 automated sequencer (Applied Biosystems). I aligned 
sequences using CodonCode Aligner v3.5.6 (CodonCode, 
Dedham, MA) and Geneious 4.8.5 (Biomatters). To ensure 
sequence accuracy, I obtained, aligned, and manually edited 
sequences from both strands (derived from each forward and 
reverse primer) belonging to the same individual. I used the 
consensus sequence of each individual for further analyses. 
I deposited consensus sequences in GenBank (JN849580- 
JN849633). 

Microsatellite genotyping 

In order to maximize the chances of detecting population 
structure across the Panama Canal, I genotyped microsatel- 
lite loci. Due to their higher mutation rates (Jarne and 
Lagoda 1996), microsatellites provide information about 
genetic structure over shorter time scales than mtDNA, 
and provide a multilocus perspective on population genetic 
structure. I amplified seven polymorphic microsatellite loci 
(Table 3) from previously published studies (Escobar- Paramo 
2000; Bohle and Zischler 2002) on ABI 2720 (Applied Biosys- 
tems) or BioRad iCycler (Bio-Rad) thermocyclers using flu- 
orescently labeled primers. Cycling conditions followed the 
mitochondrial protocol with the following modifications: 35 
cycles of amplification and locus-specific annealing temper- 
ature (see Ta, Table 3). I genotyped samples in an ABI 3730 
automated sequencer (ABI, Foster City, CA): 1 /xl of PCR 
product was added to 8.8 /xl of formamide with 0.2 /A of 
GeneScan 500-LIZ size standard (ABI). I scored genotypes 
manually using Genemapper 4.0 (ABI). To ensure robust- 
ness of genetic data, I typed homozygote genotypes from 
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Table 3. Per-population microsatellite characteristics. Population means ± SD for each are presented. 



Population 


Locus 


N 


Ta (°C) 


Na 




Ar 


Par 


Ho 




He 






UHe 


Gamboa 


Mean 


19 




2.86 ± 


1.21 


2.35 ± 0.57 


0.19 ± 0.16 


0.55 ± 


0.21 


0.50 ± 


0. 


13 


0.52 ±0.14 




SB7 


19 


54 


5 




3.0636 


0.3341 


0.421 




0.652 






0.670 




SB8 


19 


54 


2 




1.8387 


0.2213 


0.421 




0.332 






0.341 




SB19 


19 


54 


2 




1.9801 


0.0003 


0.579 




0.478 






0.491 




SB38 


19 


54 


2 




1.8763 


0.0959 


0.474 




0.361 






0.371 




Ceb10 


19 


52 


4 




3.1310 


0.2419 


0.947 




0.680 






0.698 




Ceb11 


19 


52 


3 




2.5717 


0.4081 


0.684 




0.532 






0.546 




Ceb128 


19 


52 


2 




1.9913 


0.0080 


0.316 




0.499 






0.512 


Panama W. 


Mean 


22 




4.29 ± 


1.80 


2.85 ± 0.74 


0.51 ± 0.55 


0.54 ± 


0.21 


0.60 ± 


0. 


14 


0.61 ±0.14 




SB7 


22 


54 


5 




3.2074 


0.5249 


0.455 




0.658 






0.673 




SB8 


22 


54 


7 




3.3217 


1.3316 


0.591 




0.684 






0.700 




SB19 


21 


54 


2 




1.9914 


0.0012 


0.524 




0.500 






0.512 




SB38 


22 


54 


6 




4.0426 


1 .2063 


0.955 




0.785 






0.803 




Ceb10 


22 


52 


3 




2.7673 


0.0774 


0.500 




0.624 






0.638 




Ceb11 


22 


52 


4 




2.6181 


0.2420 


0.455 




0.567 






0.580 




Ceb128 


22 


52 


3 




1.9836 


0.1592 


0.273 




0.361 






0.369 


_>uutr[ al lo 


Mf?3n 


1 8 




4 14 + 


2 34 


2 37 + o 99 


n 48 + n 58 

U.^+O _!_ U.JO 


0 61 + 


0 1 9 


0 57 + 


o 


1 6 


n QQ + n 1 7 




SB7 


17 


54 


8 




4.2180 


1.4753 


0.882 




0.782 






0.806 




SB8 


16 


54 


6 




3.6544 


1.0413 


0.750 




0.689 






0.712 




SB19 


17 


54 


2 




1.9356 


0.0001 


0.588 




0.415 






0.428 




SB38 


18 


54 


5 




3.7931 


0.4310 


0.722 




0.748 






0.770 




Ceb10 


18 


52 


4 




2.5491 


0.4001 


0.444 




0.451 






0.463 




Cebll 


18 


52 


2 




1.9584 


0.0142 


0.333 




0.444 






0.457 




Ceb128 


17 


52 


2 




1.9544 


0.0004 


0.529 




0.438 






0.451 



N = number of individuals typed; Ta = annealing temperature for reaction; Na = number of alleles loci used in this study; Ar = allelic richness. Par = 
private allelic richness; Ho = observed heterozygosity; He = expected heterozygosity; UHe = unbiased heterozygosity. 



the museum skin samples and field collected hair samples 
from at least two independent PCR reactions. I genotyped a 
subset of samples from the Gamboa population from both 
hair and tissue samples. Additionally, a subset of samples 
was genotyped de novo from independent extractions. I de- 
posited microsatellite genotypes in the Dryad data repository 
(doi: 10.506 l/dryad.mg87590q) . 

I checked loci for evidence of null alleles and geno typing er- 
rors using the program Microchecker 2.2 (van Oosterhout 
et al. 2004). I tested microsatellite loci for deviations from 
Hardy-Weinberg and linkage disequilibrium using Fstat 
2.93 (Goudet 1995). 

Genetic diversity 

I used the program DnaSP 5.1 (Librado and Rozas 2009) 
to calculate nucleotide (jr) and haplotypic (h) diversity for 
mitochondrial sequence data. I calculated number of alleles, 
expected and observed heterozygosities for microsatellite loci 
using Arlequin 3.1 and GenAlex 6.1 (Peakall and Smouse 
2006). I calculated allelic richness and private allelic rich- 
ness using the rarefaction method (to control for sample size 
differences) in Hp-Rare 1.1 (Kalinowski 2005). I tested for 
differences between populations in microsatellite genetic di- 
versity statistics using analysis of variance (AN OVA). Statisti- 



cal tests are two-tailed and I report means with their standard 
deviations (mean ± SD), unless otherwise noted. 

Population genetic structure 

I examined mitochondrial haplotype relationships using a 
parsimony network as calculated by Tcs 1.22 (Clement et al. 
2000), with a 95% connection limit. I conducted an exact test 
of population differentiation based on mitochondrial haplo- 
type frequencies (Raymond and Rousset 1995) in Arlequin 
3.1 (Excoffier et al. 2005). I calculated pairwise F-statistics 
in Arlequin to investigate population genetic structure via 
mtDNA and microsatellite data. I tested for differences in 
Fst values between populations using 10,100 permutations. I 
used analysis of molecular variance, AM OVA, (Excoffier et al. 
1992) to examine the amount of genetic variance explained 
by within- and among-population (i.e., sampling locality) 
variation. 

Population assignment 

To determine the probable number of populations repre- 
sented by the dataset, I used three different Bayesian cluster- 
ing methods, due to previously reported variability in perfor- 
mance of these methods (Waples and Gaggiotti 2006; Rowe 
and Beebee 2007) and the somewhat different algorithms 
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implemented and information used in each program. I ran 
Structure (Pritchard et al. 2000) that was run with no pre- 
vious population information (USEPOPFLAG = 0), using 
an admixture model and assuming allele frequencies were 
correlated. I conducted four replicate runs at k = 1-5 using a 
burin of 10 6 and a run length of 10 9 . 1 selected the optimal k 
value using the highest Pr(Xlk) values (Pritchard et al. 2000). 
I conducted an assignment test in Geneclass2 (Piry et al. 
2004) to determine the probability that each individual was 
assigned to its own "population," in this case representing the 
sampling locality. I also ran the analysis with an alternate pop- 
ulation definition collapsing the Soberania and Panama West 
individuals into one population. Individuals were assigned 
to a population using the criterion of Rannala and Mountain 
(1997) and the probability of these assignments was calcu- 
lated using a Monte-Carlo resampling technique (Piry et al. 
2004), based on 10,000 simulated individuals. The Type I er- 
ror was set at 0.01. I used Geneland as a third method of 
population clustering because it incorporates spatial data in 
order to identify genetic discontinuities in a spatially explicit 
fashion, which is relevant to the question at hand. First, I 
determined the most probable number of k populations us- 
ing the Markov chain Monte Carlo (MCMC) over 5 x 10 5 
iterations (Coulon et al. 2004), using the uncorrelated model. 
I conducted five replicates of this process to ensure consis- 
tency of results (Rowe and Beebee 2007). I generated spatial 
maps of the posterior probability of population membership 
by using the posterior probability of population membership 
obtained from the MCMC simulation. I overlaid the course 
of the Chagres River and the Panama Canal (extracted from a 
satellite photo) was overlaid on the population map, in order 
to investigate the coincidence of geographic barriers with the 
population limits calculated by Geneland. 

Results 

Genetic samples were collected from a total of 59 S. geoffroyi 
across the Panama Canal watershed. The number of indi- 
viduals from which microsatellite (and mitochondrial) data 
were collected in each population were: Gamboa: 19 (17), 
Panama West: 22 (21), and Soberania: 18 (16). 

Microsatellite genotyping 

Microsatellite loci showed no evidence of null alleles or geno- 
typing errors. Loci did not show evidence of departure from 
Hardy-Weinberg equilibrium or linkage disequilibrium in 
any of the populations after correction for multiple tests at 
the 0.05 nominal level. Genotyping of hair samples and ear 
tissue samples from Gamboa yielded identical genotypes at 
all loci, suggesting that genotypes from hair samples were 
reliable. 



Table 4. Diversity statistics for mitochondrial sequence data. Standard 
deviation is for sampling and stochastic processes. 







Number 


TV 




h 






of 


Nucleotide 




Haplotype 


Population 


N 


haplotypes 


diversity 




diversity 


Panama West 


21 


8 


0.00289 ± 


0.00047 


0.886 ± 0.036 


Gamboa 


17 


2 


0.00392 ± 


0.00033 


0.529 ±0.045 


Soberania 


16 


7 


0.00171 ± 


0.00026 


0.792 ± 0.076 



N = number of individuals sequenced in each population. 



Genetic diversity 

Analysis of mitochondrial sequences revealed 12 variable sites 
and 13 haplotypes. As expected for a larger geographic sam- 
ple, the Panama West population had the largest number of 
haplotypes, but it did not have the greatest nucleotide diver- 
sity, as did Gamboa (Table 4). Measures of allelic diversity 
at microsatellite loci (Table 3) were not significantly differ- 
ent between sampling localities (ANOVA's: Observed het- 
erozygosity P = 0.8008, unbiased heterozygosity P = 0.51 18, 
allelic richness P = 0.3991, rarefaction-calculated private al- 
lelic richness P = 0.3882). Gamboa appeared to have the 
lowest allelic diversity across measures, but these differences 
were not statistically significant. The lack of difference in 
microsatellite diversity included the Panama West locality, 
despite the differences in the nature of the sampling regime 
(larger geographic extent, different time frame). 

Population structure 

The haplotype network (Fig. 2) revealed sharing of haplo- 
types between the Panama West and Soberania sites. The 
number of haplotypes observed in the Panama West locality 
was larger than at the other localities. The Gamboa local- 
ity was composed of two haplotypes, which were not shared 
with the other populations. An exact test based on haplotype 
frequencies suggested strong evidence for differentiation of 
three populations (P < 0.0001; 30,000 Markov steps). 

Permutation analyses of Fst values yielded significantly 
different (P < 0.001) values between all pairwise compar- 
isons of sampling localities (Table 5). All Fst values were 

Table 5. Pairwise comparisons for F-statistics using mtDNA and mi- 
crosatellite data as calculated by ARLEQUIN. Microsatellite data are 
above the diagonal and shaded. Statistically significant differences at 
"P < 0.001 and "P < 0.0001, respectively. P values calculated based on 
10,100 permutations. 



Population 


Panama West 


Gamboa 


Soberania 


Panama West 




0.13616* 


0.06428* 


Gamboa 


0.28532* 




0.13247* 


Soberania 


0.12061** 


0.34123" 
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Table 6. Analysis of molecular variance for microsatellites and mitochondrial (in parentheses) data. Based on 10,100 permutations (P < 0.0001). 



Source of variation 


df 


Sum of squares 


Variance components 


Percentage of variation 


Among populations 


2(2) 


22.95 (24.47) 


0.24 (0.60) 


11.11 (27.85) 


Within populations 


115(51) 


223.70 (79.00) 


1.95 (1.55) 


88.89 (72.15) 


Total 


117 (53) 


246.703 (103.46) 


2.19 (2.15) 




Fixation index 


Fst: 0.11 (0.28) 









statistically significantly different from zero. Fst values cal- 
culated from mitochondrial sequence data were in general 
over two times larger than those calculated from microsatel- 
lite data. In each case, Fst values calculated from localities 
across the Chagres River were around two times more di- 
vergent relative to comparisons across the Panama Canal, 
suggesting greater differentiation across the older riverine 
barrier. 

AM OVA for mitochondrial data attributed 27.85% of vari- 
ance to among-group (sampling locality) variation, com- 
pared to 11.11% for microsatellite data (Table 6). The 
fixation index and among-population variance components 
calculated by both AMOVA's differed significantly from ran- 
dom expectation (P < 0.0001). 

Population assignment 

The three Bayesian methods yielded different estimates for 
most probable number of populations. Structure detected 
k = 2 as the most likely number of populations, with the 
Gamboa locality distinct from the combined Panama West 
and Soberania localities. The assignment plot is depicted in 
Figure 2. Geneclass2 analysis correctly assigned 78% of in- 
dividuals to their sampling localities when three populations 
were assumed. When two populations were assumed, assign- 
ment success increased to 88.1% of individuals as did the 
quality index (k = 3: 62.06% k = 2: 80.77%), which repre- 
sents the mean value of individual assignment scores (Piry 
et al. 2004). Both Structure and Geneclass2 recovered lo- 
calities across the Chagres as distinct populations, but did 
not always recover two distinct populations when comparing 
localities across the Panama Canal. Geneland on the other 
hand, clearly delineated three populations with minimal vari- 
ance in the posterior probabilities of population estimation 
over multiple runs. The location of the riverine barriers un- 
der study were largely consistent with the population limits 
delineated by Geneland (Fig. 4). 

Discussion 

The results suggest that both the Chagres River and the 
Panama Canal have contributed to population genetic struc- 
ture in S. geoffroyi inhabiting the Panama Canal watershed. 
Although the sampling regime is limited, the results pro- 



vide good, albeit preliminary, evidence of differentiation 
across two riverine barriers. Analyses of F-statistics, haplo- 
typic data, and output from Bayesian assignment algorithms 
are collectively consistent with the Chagres River playing 
a role in relatively strong population differentiation, espe- 
cially considering the small geographic distances between 
sampling localities (~6 km). All analyses except Structure 
and Geneclass2, indicate that there is detectable popu- 
lation differentiation among sampling localities across the 
Panama Canal. As expected, the level of differentiation was 
smaller across the Panama Canal than across the Chagres 
River. 

There were possible limitations imposed by the sampling 
regime for interpreting differentiation of sampling localities 
across the Panama Canal. In particular, the wider geographic 
sampling at the Panama West locality may cause additional 
allelic and haplotypic variation to be sampled. Analyses of 
microsatellite data suggest that this is not the case. Although 
the number of haplotypes is larger in the Panama West local- 
ity, it was Gamboa that had the greatest nucleotide diversity. 
However, because of these differences in sampling regime, 
the results of this study should be interpreted as preliminary 
evidence. More generally, the limited number of sampling lo- 
calities also underscores the need for caution absent broader 
geographic sampling. Previous studies investigating the role 
of riverine barriers have found discrepant results when sam- 
pling at different localities along riverine barriers (e.g., Patton 
et al. 1994). Future studies incorporating several populations 
across a larger area in the Canal watershed will allow quantita- 
tive tests of the role of riverine barriers in creating population 
structure in S. geoffroyi. 



Differentiation discrepancies according to 
marker type 

The degree of differentiation among sampling localities in- 
ferred using mitochondrial data was larger that that cal- 
culated with microsatellite genotypes. This was true across 
both the Chagres River and the Panama Canal. The AM OVA 
conducted on both genetic datasets indicated that a greater 
proportion (almost three times) of among-group variance 
was explained by mitochondrial sequence differences, as ex- 
pected for a uniparentally inherited marker. 
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been reported to perform poorly at detecting populations 
with low differentiation (Waples and Gaggiotti 2006). More- 
over, the creators of Structure caution that the large pa- 
rameter space complicates the selection of k (Pritchard et al. 
2000). This situation seems applicable to the current study 
as evidenced by the overlapping variances of Pr(Xlk) for two 
and three populations (Fig. 3). On the other hand, Geneland 
consistently identified three populations and the geographic 
projection of population membership probabilities roughly 
coincided with the approximate location of both putative bar- 
riers under study. These results underscore the variability of 
k estimates from different population clustering algorithms 
and suggest that future researchers should use multiple meth- 
ods (Rowe and Beebee 2007), and evaluate results in light of 
the biological significance to the study species ( Pritchard et al. 
2000). 



Figure 3. Posterior probability of each k estimate of Structure. Error 
bars are the variance of the posterior probability estimate. Inset is the 
bar plot for the most likely number of populations (k = 2), which shows 
the fractional assignment probability to each individual to the clusters 
inferred by Structure. 

Number of distinct populations using 
Bayesian clustering 

Although the results from Structure and Geneclass2 sug- 
gest two populations in the dataset, Bayesian algorithms have 



Riverine barriers in the Panama Canal 
watershed 

The Chagres River has been associated with genetic structure 
in at least one other species. Lampert at al. (2003) showed that 
the Chagres River formed a barrier to dispersal of tiingara 
frogs as indicated by isolation by distance patterns calcu- 
lated using microsatellite markers. Evidence that the Panama 
Canal has affected gene flow in a multitude of species is more 
abundant. Meyer et al. (2009) showed that bat populations 
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Figure 4. Map of population membership as calculated by Geneland. Pixels are colored according to the modal posterior probability of population 
membership. The approximate location of Panama Canal and the Chagres River (drawn using georeferenced satellite images) are indicated in black. 
The inset shows the density of the estimate of k (number of populations) along the Markov chain. 
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inhabiting the islands created upon the flooding of Gatun 
Lake had lower genetic diversity and higher genetic differenti- 
ation than mainland populations, according to their dispersal 
abilities. Studies of freshwater fish suggest that distinct species 
assemblages existed on either side of the Cordillera Central 
on the Atlantic and Pacific Slopes of the Isthmus (Meek and 
Hildebrand 1916; Smith et al. 2004) and now exist in the 
same communities as a consequence of the aquatic connec- 
tion provided by the Canal (Smith et al. 2004). Thus, the 
Panama Canal has affected population structure in a variety 
of taxa, increasing gene flow in aquatic species and restricting 
it in some terrestrial species. 

Demographic evidence also supports the idea of reduced 
gene flow in terrestrial species as a consequence of the cre- 
ation of the Panama Canal. Intensive studies on Barro Col- 
orado Island (BCI) have shown that multiple bird species 
have become locally extinct, most likely as a cause of the 
limited dispersal across the Canal (Robinson 1999). There 
is also demographic evidence that the canal had significant 
effects on S. geoffroyi populations: in BCI, the tamarin popu- 
lation has seen decline, as observational (Enders 1939; Eisen- 
berg & Thorington 1973) and census data (Wright et al. 
2000) suggests. While habitat conversion (from secondary 
to primary forest) has been suggested as a cause of the 
decline of tamarins on BCI, the results of this study and 
those cited above suggest that the absence of dispersal and 
gene flow could have played a part in this demographic 
change. 

Conservation implications 

The lower differentiation across the Panama Canal suggested 
by this study points to only modest structure. However, in 
the absence of migrants these populations may diverge in the 
future, as has happened more clearly across the Chagres. In 
fact, divergence may be hastened across the Panama Canal, 
due to decreasing habitat availability west of the Panama 
Canal. In contrast, the Chagres headwaters are < 25 km from 
both Soberania and Gamboa populations and well forested. 
This may ensure that over time gene flow will persist across, 
that is around, this barrier, which will likely not be the case 
for the Panama Canal. This study adds to the growing body 
of literature on the effects of recent anthropogenic barriers 
on population structure and genetic diversity (Piertney et al. 
1998; Smith et al. 2004; Epps et al. 2005). The results of this 
study highlight the utility of the Panama Canal watershed 
as an ideal testing ground for questions of population struc- 
ture. Moreover, the proximity of natural protected areas to 
two rapidly growing population centers (Rompre et al. 2008) 
provide challenges for species conservation, but ample op- 
portunities for conservation-oriented biological research on 
a number of tropical species. It is hoped that the current study 
will stimulate such research. 
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